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O ■ ABSTRACT 

C/j . Context. The evolution of magnetic fields in Ap stars during the main sequence phase is presently mostly unconstrained by observation 

, ^ , ' because of the difficulty of assigning accurate ages to known field Ap stars. 

Aims. We are carrying out a large survey of magnetic fields in cluster Ap stars with the goal of obtaining a sample of these stars with 
well-determined ages. In this paper we analyse the information available from the survey as it currently stands. 
^ , Methods. We select from the available observational sample the stars that are probably (1) cluster or association members and (2) 

^— ' magnetic Ap stars. For the stars in this subsample we determine the fundamental parameters T^s, L/Lq, and M/Mq. With these data 

and the cluster ages we assign both absolute age and fractional age (the fraction of the main sequence lifetime completed). For this 
[ purpose we have derived new bolometric corrections for Ap stars. 

. Results. Magnetic fields are present at the surfaces of Ap stars from the ZAMS to the TAMS. Statistically for the stars with M > 3Mo 

• ' the fields decline with advancing age approximately as expected from flux conservation together with increased stellar radius, or 

' perhaps even faster than this rate, on a time scale of about 3 10' yr. In contrast, lower mass stars show no compelling evidence for 

I field decrease even on a timescale of several times 10** yr. 

. Conclusions. Study of magnetic cluster stars is now a powerful tool for obtaining constraints on evolution of Ap stars through the 

■ main sequence. Enlarging the sample of known cluster magnetic stars, and obtaining more precise RMS fields, will help to clarify the 

results obtained so far. Further field observations are in progress. 
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1. Introduction all have large-scale (global) magnetic fields with strengths in 

the range of a few hundred to about 30000 G. Their distinc- 
A small fraction, of the order of 10 %, of the main sequence A ^^^^ chemical anomalies appear to vary fairly systematically with 
and B stars show obvious signs of anomalous atmospheric chem- j^^^. ^^^^^^ ^^^^ s^^rEu stars with 7000 < T.ff < 10000 
ical abundances in their spectra. The distinctive surface abun- ^^^^ ^^^^^^ overabundances of Sr, Cr, and rare earths, while 
dances of these stars have never found a satisfactory explana- ^^^^^^ ^^^^ ^^^^^ ^j^^ 10000 < < 
tion as the result of formation in a chemically anomalous region ^5000) typically show strong Si overabundance and He under- 
of the galaxy, or as the product of internal nuclear evolution or abundance. The hottest Ap stars (the He-str stars with ~ 
of surface nuclear reactions. Instead, it is now generally thought 20 000 K) have overabundant He. The broad features of which 
that these stars have undergone chemical fractionation as a result eiej^ents are expected to be over or underabundant as a func- 
of competition between gravitationally driven diffusion of trace ^^^^ 7.^^ ^.^ understood - to some extent - as a result of the 
elements, radiative levitation of some of these elements, and in- competing processes mentioned above, 
teraction of the separation process with mass loss and with mix- 
ing processes such as convection and turbulence (Michaud |1970l The magnetic Ap stars are also usually periodically variable, 
Babel [T992l Richer et al.l2 000l ). in light, spectrum, and magnetic field. The periods range from 
One of the major groups of these main sequence stars with about half a day to decades without any obvious difference in the 
anomalous atmospheric abundances is the "magnetic peculiar A" phenomena of variation, and the periods observed are inversely 
stars, which are commonly refeiTed to as Ap stars. These stars correlated with Vg sin /. It is clear from these facts that the period 

of variation must be the rotation period. This variability implies 

Send offprint requests to: J. Landstreet that the field and the chemical abundances are non-uniformly 



2 



J. D. Landstreet, S. Bagnulo, V. Andretta, et al.: The evolution of magnetic fields in Ap stars 



distributed over the surface of an Ap star. The rotation axis is 
not an axis of symmetry, so that the observer sees a configuration 
that changes as the star rotates. 

Much is now understood about these magnetic Ap stars. 
Simple models of field structure indicate that almost all mag- 
netic Ap stars have roughly dipolar field structures. The axis of 
(approximate) symmetry of the field is usually oblique to the 
rotation axis, so that the line-of-sight component of the field 
varies roughly sinusoidally with time (e.g. Mestel & Landstreet 
I2005I I. Often several chemical elements are distributed quite 
non-uniformly over the surface, and low-resolution models (e.g. 
Strasser, Landstreet & Mathys 2001) of the distributions are 
available for a number of Ap stars. It is now becoming possi- 
ble to deduce quite detailed maps of both field structure and el- 
ement distribution for a few stars for which the variations have 
been observed spectroscopically in all four Stokes parameters 
(Kochukhov, Bagnulo, Wade et al. [20041 KochukhovjIOOl. 

The lack of any energetically significant convection zone 
near the surface of all but the coolest Ap stars, and the absence 
of a direct correlation between rotation rate and field strength (as 
is found in solar-type stars, see Mestel & Landstreet 2005 ) has 
led to general agreement that the magnetic fields observed are 
probably fossil fields retained from an earlier evolution stage, a 
situation which is possible because of the extremely long decay 
times. 

Numerical calculations have suggested how a fossil mag- 
netic field might evolve with time during the long main sequence 
phase as a consequence of ohmic decay, changes in stellar struc- 
ture, and internal circulation currents (Moss l200TT l. but the rela- 
tionship of such computations to observed stars is still very un- 
clear. Braithwaite & Spruit (2004) and Braithwaite & Nordlund 
(I2OO6I 1 have shown that the combination of the Cowling insta- 
bility with rotation and meridional circulation may lead to major 
reorganisation of the internal flux. Similarly, theoretical com- 
putations suggest how atmospheric chemical peculiarities might 
develop and evolve with time through the main sequence phase, 
but so far these computations have not succeeded in reproduc- 
ing observed abundance patterns in detail (Babel & Michaud 
[1991 ). Observational data on fields and abundances in a sample 
of Ap stars for which the absolute age and fraction of the main 
sequence evolution completed (we will call this the "fractional 
age") is known would be extremely valuable to test theoretical 
ideas and provide clues to important processes. 

At present, we have few constraints from observations about 
the evolution of either fields or atmospheric chemistry with time 
before or through the main sequence phase. Ap stars have some- 
what anomalous, and variable, energy distributions, which make 
it difficult to place these stars in the HR diagram precisely 
enough to obtain accurate ages for individual stars from com- 
parison with theoretical evolution tracks. The problem is com- 
pounded by the fact that the atmospheric chemistry does not 
readily reveal the bulk stellar abundances, and so we do not even 
know how to choose the metallicity used in computing the evo- 
lution tracks which are to be compared with observations. As a 
result, until recently very few Ap star ages could be securely de- 
termined. Age determinations of field stars are barely accurate 
enough to determine whether a given star is in the first half or 
second half of its main sequence life. (This problem is discussed 
in more detail by Bagnulo et al. 120061 hereafter Paper I). 

The problem of obtaining an accurate age for an Ap star is 
considerably simpler if the star is a member of a cluster. Since 
the age of a cluster can be determined typically to somewhat 
better than +0.2 dex (by isochrone fitting to the main sequence 
turnoff or to low-mass stars evolving to the zero-age main se- 



quence, or by studies of Li depletion in low-mass stars), we can 
use the cluster age to determine the fraction of the main sequence 
lifetime which has elapsed. For a star in a cluster whose age is 
small compared to the star's main sequence lifetime, this evolu- 
tionary age is relatively precise. On the other hand, age resolu- 
tion is not very good for the second half of a star's main sequence 
lifetime (see Paper I, and Sect. 5.1). 

It has recently become practical to obtain magnetic measure- 
ments of a statistically interesting number of stars in clusters. In 
Paper I we have reported the first large scale observational sur- 
vey of magnetic fields in peculiar (and normal) A and B stars in 
open clusters and associations. This survey has been motivated 
by our desire to obtain empirical information about the evolution 
of magnetic fields and atmospheric chemistry in middle main 
sequence stars as these stars evolve from the zero-age main se- 
quence (ZAMS) to the terminal-age main sequence (TAMS). 

In the present paper, we draw some first conclusions from 
this (ongoing) survey about the evolution of magnetic fields from 
ZAMS to TAMS. In Sect. |2] we will examine and select all the 
stars that are probable cluster or association members and proba- 
ble Ap stars, and for which magnetic measurements are available 
(either from our own survey or previous works). In Sect. |3]we 
determine the fundamental parameters Teff and L/L© for these 
stars. Combining this information with suitable evolution tracks, 
we derive stellar masses and fractional ages for the stars of the 
sample in Sect. |4] In Sect. |5] we comment on the characteristics 
of our sample (distributions of age, mass, and magnetic field). 
Finally, in Sect. |6] we consider what these data reveal about the 
evolution of surface fields from the ZAMS to the TAMS. Sect 7 
summarises our conclusions. 

2. Cluster membership of stars with magnetic data 

In Paper I we presented new magnetic measurements of many 
candidate or confirmed Ap stars that may be members of clus- 
ters or associations (as well as a large number of measurements 
of probable cluster members that are not Ap stars). From these 
data and from field measurements in the literature, we have as- 
sembled a database of all probable Ap cluster/association mem- 
bers for which usefully precise magnetic field measurements are 
available. This database is presented in Table |2j the content will 
be explained below. In this section we analyze more critically 
the evidence for cluster membership, and for magnetic Ap na- 
ture, for the stars in Table |2] This will enable us to remove any 
stars for which either membership or Ap nature is doubtful from 
the stellar sample we finally analyse. 

2. 1 . Assessing cluster membership 

Cluster (or association) membership may be tested in several 
ways. Traditionally, membership is tested by plotting a colour- 
magnitude diagram for the cluster. Stars which lie too far above 
or below the cluster main sequence are rejected as members. 
Previous studies (e.g. North , 1993. Maitzen, Schneider & Weiss 
119881 ) have shown that Ap stars that are probable cluster mem- 
bers lie reasonably close (in observational HR diagrams, for ex- 
ample Johnson V magnitude as a function of colour B - V) to 
the main sequence of normal stars, although Ap stars sometimes 
deviate from the normal main sequence by somewhat more than 
normal cluster members. In general, the stars in Table|2]are suffi- 
ciently close to the main sequences of the clusters to which they 
may belong that this criterion does not immediately eliminate 
any of the stars in the table from membership. The HR diagram 
test of membership will be applied again later in this paper when 



J. D. Landstreet, S. Bagnulo, V. Andretta, et al.: The evolution of magnetic fields in Ap stars 



3 



the provisionally accepted cluster members are placed in the the- 
oretical {L/Lq vs Tgff) HR diagrams of their specific clusters. 

Comparisons of measured parallaxes and proper motions of 
candidate stars with the average values for the cluster provide 
powerful tests of membership. The situation for data on both 
these properties has improved dramatically as a result of the 
Hipparcos mission (ESA 1997). This mission led to measure- 
ments of parallaxes and proper motions of about 120000 stars, 
with typical errors of the order of 1 mas and 1 mas yr"', so that 
accurate (+25 %) parallaxes are available for tens of thousands 
of stars out to a distance of about 300 pc. Observations of a much 
larger number of stars than that included in the Hipparcos par- 
allax programme, using the star mapper instrument on the satel- 
hte, led to the Tycho (or Tycho-1) Catalogue (H0g et al. 1 19971 1 
of proper motions of about 1 million stars, with accuracies of 
the order of 10 mas yr"' (at V = 9) in both coordinates. The 
Tycho data were combined with positions from the Astrographic 
Catalogue to produce the Tycho Reference Catalogue (H0g et 
al. 1 19981 1. with proper motions of 990 000 stars accurate to about 
2.5 mas yr"', and an extension of this project led to the Tycho-2 
catalogue (H0g et al. 2000a 2000b ), which combines data from 
the Hipparcos mission with those from many ground-based po- 
sition catalogues to provide proper motions for about 2.5 million 
stars with typical proper motion standard errors of 2.5 mas yr 

These data have been used by several groups for im- 
portant studies of membership in clusters and OB associa- 
tions. Robichon et al. (1999 ) determined mean distances, mean 
proper motions, and membership of Hipparcos parallax stars 
located near the centres of 18 clusters less than 500 pc away. 
Membership was derived from assessment of position, parallax, 
proper motion, and photometry. The work led to list of accepted 
members of each cluster; the absence of a Hipparcos parallax 
star from one of these lists presumably means that it was not 
accepted by Robichon et al. dl 9991 1 as a member of the corre- 
sponding cluster. 

Baumgardt et al. (I2000I I derived mean distances and mo- 
tions for 205 open clusters using Hipparcos parallax stars found 
in these clusters. Membership of Hipparcos parallax stars in 
these clusters was assessed using photometry, radial velocity, 
position, and proper motion data from the ground and from 
Hipparcos, together with some additional proper motions from 
the Tycho Reference Catalogue. Baumgardt et al. report both 
an overall membership assessment (member, possible member, 
non-member) and a proper motion membership probability. 

The Hipparcos data were used by de Zeeuw et al. ( 1 19991 ) 
to assess membership of Hipparcos parallax stars in nearby OB 
associations. Several of the associations included in their study 
(Sco OB2, a Per, CoUinder 121, Trumpler 10, etc) are repre- 
sented in our survey. 

The release of the Tycho-2 proper motion catalogue made 
possible further large-scale studies of membership in open clus- 
ters. Dias, Lepine and Alessi (2001] 120021 ) have determined 
mean proper motions of more than 200 clusters, many of which 
were not included in the studies discussed above, using in most 
cases several tens of stars per cluster. They have derived mem- 
bership probabilities based only on the proper motions, which 
are tightly clustered for an open cluster but (usually) have a con- 
siderably larger dispersion for field stars. 

Kharchenko et al. (2005) created a database of possible stel- 
lar members of some 520 clusters. They use a combination of 
Tycho-2, Hipparcos, and ground-based data to determine av- 
erage cluster parameters (distance, reddening, proper motions, 
etc). Memberships of individual stars in their database are as- 
sessed using all available information (position, photometry. 



proper motions, etc), and then cluster ages are derived using the 
positions of the brightest stars relative to theoretical isochrones. 

We have used the studies discussed above together with the 
actual astrometric data to assess the probability that the stars 
of Table |2] are actually cluster members. Wherever possible we 
have adopted the consensus of the previous membership studies. 
However, the individual determinations of membership proba- 
bility differ somewhat from one study to another, both because a 
study may be based either on Hipparcos or Tycho-2 proper mo- 
tions, and because the stars selected as definite cluster members, 
from which the cluster mean motions are deduced, varies from 
one study to another. Note in particular that Kharchenko et al. 
(i2005) often assign membership probabilities based on proper 
motions that are significantly smaller than those found in any 
of the other studies (e.g. for HD 21699 in the a Per cluster, 
de Zeeuw et al. [19991 give P = 0.84 and Robichon et al. [19991 
consider the star to be a definite member, but Kharchenko et al. 
(2005) give = 0.25 as the proper motion membership prob- 
ability). We generally consider that the Kharchenko et al. data 
support membership if their proper motion membership proba- 
bility is larger than about 0. 15 (their 2-cr limit). 

In a few cases the literature conclusions on membership are 
discordant, or no membership probability based on recent as- 
trometry is available. In these cases we have compared proper 
motions and parallaxes, when available, to the cluster means, and 
assessed membership ourselves. We consider a star to be a prob- 
able member if both proper motions (and parallax if available) 
are less than 2.5cr from the mean cluster values, and a definite 
member if all are within Icr of the cluster means. If no accu- 
rate astrometry is available, but the star is within the bounds of 
the cluster and has colour and magnitude consistent with mem- 
bership, we have generally assumed that the star is a probable 
member 

In the particular case of stars in Ori OBla, OBlb, and OBlc, 
the members of the association are not well separated in proper 
motion from the field. We have adopted the criterion of Brown 
et al. (119991 ) for proper motion membership, and have tested ap- 
parent brightness and parallax for consistency with membership 
as far as possible. However, membership of stars in these regions 
of Ori OBI is uncertain (see discussion in de Zeeuw et al. 1999). 
In Ori OB 1 d, where background obscuration virtually eliminates 
background stars, we can be more confident of membership. 

The results of this exercise are reported in Table|2]in the four 
columns under the general heading "cluster membership". In 
these columns we report our conclusions as to membership prob- 
ability considering only parallax ("tt"), proper motions ("//"), 
location in the HR diagram relative to the cluster isochrone 
("phot", see below), and overall ("member?"), using a four-level 
scale of "y" = member, "p" - probable member, "?" = probably 
non-member, and "n" = non-member The "member?" column 
also contains references to earlier membership studies. Notes at 
the bottom of Table|2] signaled by "*" in the "member?" column, 
discuss particular problems with deciding on cluster or associa- 
tion membership. 

2.2. Assessing spectral type and chemical peculiarity 

We next turn to methods to determine whether a particular star 
is an Ap star of the type in which a magnetic field is almost al- 
ways detected when sufficiently precise field measurements are 
obtained (Auriere et al. , 20041 ). The reason that this is important 
is that the field measurement errors achieved in general in Paper 
I are only small enough to detect fields in about half of the mag- 
netic Ap stars we have observed, and since we will want to ex- 
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amine the distribution of field strengths in magnetic Ap stars as 
a function of age and mass, it is important to consider both the 
Ap stars in which we detect a field and those in which no field 
has yet been found. 

Several criteria are available that can indicate Ap type spec- 
tra. Spectral classification can provide an important test of Ap 
characteristics. For many of the stars in our sample, one or 
more classification spectra are available. Classification surveys 
of a number of clusters with the goal of identifying Ap stars 
have been reported by Hartoog (119761 1 19771 1 and Abt d 19791 1 
and references therein. Spectroscopic studies of single clus- 
ters have also yielded many spectroscopic classifications (these 
may be found in the WEBDA database of open cluster data 
at [http : //w w w.univie .a c . at/webda/ , described by Mermilliod & 
Paunzen 120031 on the SIMBAD database of stellar data at 
http://cdsweb.u-strasbg.fr/; and in the compilation of data on Ap 
stars in clusters by Renson 1992 ). Unfortunately, since the clas- 
sification criteria for Ap spectra are not always unambiguous, 
and may not be clearly visible in low dispersion, low signal-to- 
noise classifications spectra, the available classifications of a star 
are frequently contradictory as to whether the star is actually an 
Ap. 

Another valuable criterion for Ap type is furnished by the 
Afl photometry system developed by Maitzen and collaborators 
(e.g., Maitzen 1993 ), and by the Z index which can be computed 
for stars for which Geneva photometry is available (e. g. Cramer 
& Maeder |1979l|1980l Hauck & North|T993). Both of these pho- 
tometric indices are sensitive to the broad energy depression at 
about 5200 A which has been found to be a robust indicator of 
Ap spectral type in roughly the temperature range 8 000 - 14 000 
K (e.g., Kupka, Paunzen, & Maitzen 2003 ). Both these photom- 
etry systems provide a large amount of information about cluster 
Ap stars; Maitzen's group has surveyed a large number of clus- 
ters for Ap stars, and Geneva data exist for many cluster stars as 
well (see http://obswww.unige.ch/gcpd/). 

A third useful criterion of Ap type is the detection of periodic 
photometric variability of the type frequently found in Ap stars. 
The catalogues of Renson & Catalano (2001') provide very use- 
ful summaries of the many particular studies found in the litera- 
ture. We have searched the Hipparcos catalogue (ESA |1997] l for 
new photometric variables, although this provided only a very 
small number of variables that had not already been identified 
by ground-based photometry. 

Finally, the definite detection of a magnetic field is a clear 
indication of Ap classification. 

We have summarized our view as to whether the star is an 
Ap star in the four columns under the general heading "Ap star" 
of Table |2] We have used the same scale adopted to evaluate 
the membership ("y", "p", "?", and "n"; a blank space denotes 
absence of information). In the column labelled "Sp" we report 
the conclusions obtained from the available spectral classifica- 
tions. We have adopted the criterion that Aa photometry defi- 
nitely supports Ap classification if Aa > 0.014, and that Geneva 
photometry supports Ap classification if Z < -0.016 (within the 
temperature range where photometric detection of peculiarity is 
valid). The conclusions are given in column Aa , Z. The degree 
of support for Ap classification from photometric variability is 
summarized in the column labelled "var". The column "mag fid" 
reports our assessment as to whether the star has a detected field, 
and also provides references to field measurements. From this 
column it will be clear that a large fraction of all field measure- 
ments of open cluster stars have been carried out in the course of 
our survey (Paper I; measurements denoted by "F" in the "mag 
fid" column); the only substantial previous surveys have been 



restricted to the Sco OB2 and Ori OBI associations. Finally, an 
overall assessment as to whether the balance of evidence sup- 
ports an Ap classification is provided in the last column "Ap?" 
of Table |2] Note that a significant number of stars are given a 
"y" or "p" even though no magnetic field has yet been detected 
in them. 



3. Determining the physical properties of cluster Ap 
stars 

With the analysis described in the previous Sect, we have iden- 
tified a subsample of about 90 stars that are definite or probable 
("y" or "p") cluster/association members, definite or probable 
Ap stars, and for which at least one magnetic field measurement 
is available (regardless of whether a field has been detected or 
not). 

For each star of this subsample we need to determine age 
and mass. These will be calculated in Sect.lH using suitable evo- 
lutionary models, assuming that we know the cluster age, the 
star's temperature and the star's luminosity. In this Sect, we de- 
scribe how we determine these physical properties. Furthermore, 
to each star we associate a value representative of the star's field 
strength. 

3.1. Cluster ages 

Because our essential goal is to provide a sample of magnetic 
measurements of Ap stars of reasonably well known ages, the 
ages of the clusters to which the stars in our sample belong must 
be critically discussed. Cluster ages for a large sample of clus- 
ters, including essentially all those considered here, are available 
from the WEBDA database (2003). Ages in this database are de- 
rived from a survey of literature values carried out by Loktin et 
al. (2001 ). These ages are in most cases based on isochrone fits 
to dereddened cluster main sequences. 

Another extensive study of ages has been carried out by 
Kharchenko et al ("2005 ) based on databases of potential clus- 
ter members generated from large observational resources such 
as the Hipparcos and Tycho databases. The ages of Kharchenko 
et al are derived by fitting isochrones individually to the most 
massive main sequence stars found among their own selection 
of cluster members, rather than by isochrone fitting to the up- 
per end of the cluster main sequence as a whole. The ages of 
Kharchenko et al. are based in many cases on only one star per 
cluster Nevertheless, these ages are in good overall agreement 
with those of Loktin et al. (1200 11 1: Kharchenko et al. have shown 
that (excluding a small number of clusters where significant dif- 
ferences in choices of cluster members lead to substantially dif- 
ferent ages) the ages from these two sources exhibit an RMS 
difference in cluster age (log f) of about cr = 0.20 dex. Note that 
neither WEBDA nor Kharchenko et al. provide uncertainties for 
the assigned cluster ages; one of our main tasks is to estimate 
these uncertainties. 

We have included these ages in our own database of possible 
cluster Ap stars, together with many recent determinations from 
the literature. After examining the available recent age determi- 
nations, we have selected what seems to us to be an appropriate 
age for each cluster, and assigned an uncertainty to this value. 
In assigning ages, we have considered both accuracy claimed by 
individual age studies and the overall agreement of various stud- 
ies. For a number of the clusters of interest, several age deter- 
minations are available, and in many cases these are reasonably 
concordant. For such clusters we have assigned an uncertainty 
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which in the best cases is taken to be ±0.1 dex. For clusters with 
fewer age determinations, or more discordant ones, we have gen- 
erally assigned an age and an uncertainty which approximately 
cover the range of the best age determinations within the range 
0.1-0.2 dex. For clusters for which no ages seem to be available, 
we have mostly used the ages from WEBDA (those of Loktin et 
al. 1200 IK and assumed an uncertainty of +0.2 dex, in agreement 
with the overall concordance between their results and those of 
Kharchenko et al. ( l2005b . 

In addition, several of the probable Ap stars we have ob- 
served turn out to be near the TAMS in the HR diagram. By 
requiring the cluster isochrone to intersect the error oval of such 
a relatively evolved star, in some cases we are able to constrain 
the cluster age more precisely than by using ages obtained from 
the literature; we obtain uncertainties that may be as small as 
0.05 dex. In these cases we have adopted the cluster ages and 
uncertainties from our isochrone fits. The procedure adopted is 
discussed further in Sect. 4. 

3.2. Stellar effective temperatures 

3.2.1. Photometry 

For stars for which Geneva photometry is available, we have 
used the FORTRAN code described by Kiinzli, North, Kurucz & 
Nicolet ( 1997 ), which is available from CDS, to determine a first 
effective temperature based on the calibration for normal (metal- 
licity 0) stars. This procedure requires supplying a colour excess 
E(B2 - Vi). We have obtained this quantity from the Johnson 
colour excess E{B - V) (either for the cluster as a whole, as re- 
ported by WEBDA or a particular study, or for the star individu- 
ally after de-reddening in the Johnson U-B,B-V colour-colour 
diagram), using the relationship E{B2 - Vi) ~ 0.88£'(B - V) 
(Hauck & North [1993). When the colours of a star are near the 
boundary between different methods of temperature determina- 
tion (e.g. XY or pTpG) we have computed the temperature us- 
ing both methods. The effective temperature determined using 
normal star calibrations is then corrected to the Ap stars tem- 
perature scale as described by Hauck & Kiinzli (119961 1. For stars 
classified as He-weak or He-strong, we have followed the sug- 
gestion of Hauck & North to adopt the temperature produced by 
the Kiinzli et al FORTRAN programme without correction. 

For stars for which Stromgren uvby/3 photometry is 
available, we have computed an effective temperature us- 
ing the FORTRAN code "UVB YBETANEW" of Napiwotzki, 
Schonberner & Wenske (|1993l l, which is in turn based mainly 
on the calibrations of Moon & Dworetsky ( 119851) . When the star 
has colours near the boundary of different calibrations (for ex- 
ample near AO), we have computed the temperature using both 
calibrations. These raw temperatures have been corrected for Ap 
stars as described by Stfpien & Dominiczak (1989) . Again, we 
have not corrected temperatures of He-weak or He-strong stars 
from the values produced by the programme. 

If only Johnson U BV photometry is available, we have de- 
reddened the star to the main sequence line in the U - B, B - V 
colour-colour diagram, and then corrected the value following 
the same prescription Stgpieri & Dominiczak ( 119891 ) use for 
Stromgren photometry. 

3.2.2. Adopted temperatures and uncertainties 

For most of the stars in our sample, both uvby/3 and Geneva pho- 
tometry are available, and we have adopted a suitable average 
of the computed temperatures. We find that temperature deter- 



mination using the two different kinds of photometry, and some- 
times more than one method of temperature determination, have 
a scatter of the order of 2-300 K. Since both methods have been 
calibrated using mostly the same "fundamental" stars, this is cer- 
tainly a lower bound to the actual uncertainty of the results. 

The global temperature calibration for Ap stars is still not 
very satisfactory. No Ap stars have a fundamental calibration 
(based on integrated flux and angular diameter measurements); 
all calibrations are based on corrections to the normal star cal- 
ibrations using a small sample of stars for which effective tem- 
peratures have been determined using the infrared flux method 
(e.g. Lanz 1985; Megessier 1988J Moni er lT992 ). fitting energy 
distributions (e.g. Adelman et al. 119951 ), or the reasonable but 
uncertain correction method used by St§pien & Dominiczak 
(i1989D . who estimate the uncertainty of their individual values 
of Teff to be 6-700 K. Furthermore, the accuracy of these cor- 
rections has recently been called into question by the results 
of Khan & Shulyak (2006j), who have calculated model atmo- 
spheres with enhanced metal content (so far generally scaled so- 
lar abundances) and the effects of a magnetic field. They find 
that, even with a rather large field and lOx solar abundances, the 
Paschen continuum slope is very similar to that of a solar com- 
position, non-magnetic star, and that the calibrations for normal 
stars recover approximately the assumed effective temperature. 

Finally, there are certainly substantial star-to-star variations 
in atmospheric composition, which make any calibration suit- 
able only on average, and not exact for any particular star This is 
clearly seen in the scatter shown in Figs. 9 and 10 of Napiwotzki 
et al. (119931) . where individual calibration points are shown along 
with the mean calibration curves adopted. 

Considering the internal uncertainties in determining the ef- 
fective temperature of a magnetic Ap star, and the further global 
uncertainties, we adopt a uniform uncertainty of +0.02 dex 
(roughly ±500 K) for our temperatures, and even this value may 
be somewhat optimistic. 

3.3. Stellar luminosities 
3.3.1. Cluster distances 

In the case of the nearest clusters (out to about 300 pc) accurate 
distances are best derived using Hipparcos parallaxes averaged 
over known members. Distances obtained in this way generally 
have uncertainties of about 0.2 mag or less in both absolute and 
apparent distance modulus (cf Robichon et al. 119991 ). 

For more distant clusters, the most accurate distances are 
obtained by isochrone fitting to dereddened cluster main se- 
quences. For a number of nearby clusters, it is possible to com- 
pare distances obtained from isochrone fitting with those ob- 
tained from parallaxes (Robichon et al l 19991 ). This comparison, 
and comparison of different determinations using isochrone fit- 
ting, allow us to estimate that a typical uncertainty for the more 
distant clusters is roughly 0.2 mag in distance modulus, similar 
to the value for the nearer clusters. In almost all cases the uncer- 
tainty in the distance modulus is dominated by the uncertainty in 
the distance rather than that in the reddening. 

We have adopted the cluster or association distances for all 
of the stars in our sample, rather than individual distances even 
when these are available from parallax observations; on average 
the mean cluster distances from parallaxes have standard errors 
that are smaller by a factor of two or more than those of indi- 
vidual cluster stars. Of course, for stars in more distant clusters, 
the best available distance is the photometric distance modulus 
of the cluster as a whole. 
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Fig. 1. Bolometric corrections using integrated fluxes of mag- 
netic Ap stars from Lanz ( 119841 red circles), Stgpeiil & 
Dominiczak (1989 blue circles), and Monier (119921 blue circles) 
as functions of effective temperature. The upper solid line is the 
BC for normal stars from Schmidt-Kaler (Lang 1992); the mid- 
dle line is Balona's (1994) fit to the fundamental data of Code et 
al. d 19761 1 and Malagnini et al. (119861 ). and the lowest curve is a 
polynomial fit to the Ap data as described in the text. 



3.3.2. Bolometric corrections 

The results of Lanz ( |1984| l strongly suggest that the bolomet- 
ric corrections for magnetic Ap stars are different from those 
of normal stars. Since we have eff'ective temperatures obtained 
from several photometric systems, we need bolometric correc- 
tions suitable for Ap stars as a function of effective temperature, 
rather than as a function of a specific photometric index (as usu- 
ally provided). We have re-examined the available data to obtain 
accurate bolometric corrections in the form we need. 

For this purpose, we have collected all the magnetic Ap stars 
for which integrated fluxes are available. Bolometric corrections 
for magnetic Ap stars have been reported by Lanz ( |1984) . while 
Stf pien & Dominiczak d 19891 1. and Monier ( |1992| l have provided 
further integrated fluxes. Almost all the stars discussed by Lanz 
are rather hot (Si and He-pec stars), and most have distances 
between 100 and 250 pc, so the observed fluxes are significantly 
affected by interstellar absorption. Except for three cool Ap stars 
in his sample, the fluxes reported by Lanz have been corrected 
approximately for extinction. The integrated fluxes reported by 
Stf pien & Dominiczak do not seem to be corrected for absorp- 
tion, even though this sample includes some (hot) Ap stars as far 
away as 150 - 180 pc, for which interstellar absorption is prob- 
ably significant. (However, for the five stars in common with 
Lanz's sample, with distances in the range 90 - 180 pc, the av- 
erage difference in integrated flux between the two studies is 
only about 4 %.) The three values of integrated flux reported by 
Monier are not corrected for interstellar absorption, but all refer 
to stars located between 50 and 100 pc from the Sun, and any 
corrections should be at most a few percent. 

We have assigned effective temperatures to the stars for 
which integrated fluxes are available using the same methods 
discussed above, and computed the bolometric corrections for 
the stars of Stfpien & Dominicz ak (19 89), and Monier (1992) 
using the first equation in Lanz ( |19841 l. Because of the lack of 
correction for interstellar extinction, we have discarded all the 



stars from Stgpieil and Dominiczak's sample that have T^ff > 
1 1 000 K. However, we verified that including them did not sig- 
nificantly change the results. 

Figure [T] shows the resulting data set of bolometric correc- 
tions as a function of T^ff together with our best fit (lower curve) 
to the Ap star mean correction BCAp(7"eff): 

BCAp(reff) = -4.891 + 15.147 61- 11.5176'^ (1) 

where = 5040.0/reff. 

We also plot the variation of bolometric correction BC„{Tsff) 
with Tf.ff for normal stars, using (upper curve) Schmidt-Kaler's 
values (Lang 1992), and Balona's 1994 fit (middle curve) to the 
fundamental data of Code et al. (119761 ) and of Malagnini et al. 
([1986 ^. 

Equation ([TJ describes the available bolometric corrections 
for Ap stars between about 7500 and 18000 K. The RMS devi- 
ation of the individual Ap star data around this curve is about 
0.1 mag, considerably larger than that for normal stars around 
Balona's fit. This reflects both larger uncertainty in both coor- 
dinates, and probably also real and substantial star-to-star varia- 
tions in the bolometric correction for these rather diverse objects. 

It is clear that the values of Code et al. (1976) and Malagnini 
et al. ( 119861 ) are the most accurate data for the normal star bolo- 
metric correction BCniTgff). Note that the BC values of Schmidt- 
Kaler, which are sometimes used for Ap stars (e.g., Hubrig et al. 
[2000) are systematically too negative compared to the more ac- 
curate data of Code et al. ( 1976) and Malagnini et al. ( 119861 ). and 
differ from the calibration adopted here by about 0.25 mag. 

With typical distance modulus uncertainties of 0.2 mag, an 
uncertainty in the bolometric correction of about 0. 1 mag, and 
reddening uncertainties of less than 0. 1 mag, we estimate that 
the typical uncertainty in Mboi is about 0.25 mag, corresponding 
to an uncertainty in log(L/LQ) of about 0.1 dex. We adopt this 
value as the standard error of luminosity values in Table [3] 

The apparent V magnitudes and values of T^ff of stars in our 
subsample have been used, together with the bolometric correc- 
tions of Eq ([1]), to determine the stellar luminosities L/Lq. For 
most clusters we used the mean cluster reddening as tabulated by 
WEBDA or one of the specific references. However, we found 
that using the mean cluster or association reddening led in sev- 
eral cases to stars appearing seriously underluminous compared 
to the position of the expected main sequence. In these cases 
we have determined the reddening E{B - V) directly for each 
star from its position in a colour-colour diagram, and used these 
values (with Ay = 3.1E(B - V)) to determine the absolute mag- 
nitude of each individual star. This was necessary for Ori OBI, 
the Upper Sco region of Sco OB2, and for NGC 2244, 2516, 
and 6193. In general, this led to a much tighter observed main 
sequence. 

3.4. Magnetic field 

Magnetic fields of Ap stars may be determined through the mea- 
surement of the Zeeman splitting of spectral lines in a simple 
intensity spectrum, or through the analysis of circular polariza- 
tion of spectral lines. The former method, exploited by Mathys 
et al. ( 119971 ). gives the average of the magnetic field modulus 
over the visible stellar disk. Zeeman splitting may be detected 
only under quite special circumstances, i.e., Ve sin / at most a few 
km s ', and field strength at least 2 kG, that are not met in most 
of the known magnetic Ap stars. Circular polarization measure- 
ments are generally a far more sensitive and broadly-applicable 
method of field detection, and give the so called mean longitudi- 
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nal magnetic field {B^), i.e., the component of the magnetic field 
along the line of sight averaged over the stellar disk. 

Because of geometrical reasons, both mean field modulus 
and mean longitudinal field change as the star rotates. In partic- 
ular, a (B,) measurement could be consistent with zero even in a 
star that possesses a strong magnetic field. For the present work 
we need to associate to each star a field value that is representa- 
tive of the stellar field strength. 

For the few stars in our sample for which the full variation 
of the longitudinal field (B ) is known and has been fitted with a 
sine wave of the form 

BAf) = Bo + Bi sin(wf - 0o), (2) 

we give the root-mean square longitudinal field strength Buns - 
(Bg + B^/iy^^. For the great majority of the stars in the table, 
however, only one measurement, or at most a few scattered mea- 
surements, are available. In this case we simply report the raw 
value of Brms averaged either over all available field measure- 
ments, or over the subset with the smallest uncertainties. Since 
most of the stars in this table have only one or two measure- 
ments, this is clearly not a very accurate statistic for an individ- 
ual star, but averaged over our sample we believe that it akeady 
provides a significant amount of useful information. 

When the individual field strength values used in the compu- 
tation of Bnns are taken from Paper I, the value derived from the 
full spectrum is used if the measurement detects a field of less 
than 1 kG. However, as discussed in Paper I, the metal lines in 
FORS 1 spectra yield an underestimate of the actual field strength 
for larger fields. Accordingly, for fields above 1 kG we have used 
only the field strength deduced from the Balmer lines. 

Note that values of Brms smaller than about 250 G typically 
should be read as indicating that the RMS field estimated from 
the current data set is of the order of 200 G or less, regardless of 
the actual tabulated value. 



4. Determining stellar masses and fractional ages 

With efifective temperatures and luminosities for almost all the 
stars in Table |3] we can now make comparison with evolution 
tracks to determine approximate masses and fractional ages (re- 
call that the fractional age of a star is the fraction of its main 
sequence lifetime, measured from the ZAMS to the TAMS, that 
has already elapsed). 

As discussed in Paper I, there is a spread of chemical com- 
positions among nearby open clusters and associations of about 
a factor of 2.5 in metallicity, so that the accurate placement of 
a star on an evolution track requires either the use of evolution 
tracks appropriate to the particular star or stellar group, or con- 
sideration of the additional uncertainties introduced by variation 
of evolution tracks with metallicity. For studies of field stars, 
for which ages are determined from comparison of position in 
the HR diagram with evolution tracks, we showed in Paper I 
that very important age uncertainties arise from this eff'ect, espe- 
cially for stars in the first half of their main sequence lives. We 
also showed that it is only currently possible, with present un- 
certainties in Teff and L/Lq, to discriminate between field stars 
in the first and second halves of their main sequence lifetimes. 
In contrast, uncertainty in bulk chemical composition does not 
introduce an important additional uncertainty into the determi- 
nation of the stellar mass. 

In this study we are determining stellar ages from the ages 
of the clusters to which the stars of our sample belong. The prin- 
cipal uncertainty for most of the available literature ages (0. 1 - 
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Fig. 2. Comparison of computed positions of stars in the Lower 
Centaurus-Lupus, Upper Centaurus-Crux, and Upper Sco re- 
gions of the Sco OB2 association with evolution tracks for stars 
of Z = 0.02 and M = 2,2.5,3,4,5,6,7, and 9 Mq, and with 
bracketting isochrones for log ages of 6.9 and 7.2, appropriate 
for these regions of the association. 

0.2 dex in logf, as discussed above) arises from uncertainty as 
to where to place the cluster turnoff" in the HR diagram, or in 
the appropriate models for pre-main sequence stars, rather than 
from uncertainty in the composition that should be adopted for 
the evolution tracks used for age determination. Consequently, 
we do not need to use composition-specific evolution tracks for 
mass and fractional age determination; we simply include the 
small uncertainty carried by this effect in our overall uncertainty 
estimates. 

We then proceed as follows. Each star is placed in the HR 
diagram, using the values of Teff and L/Lq determined previ- 
ously. Around each star we construct an approximate error oval, 
using the adopted uncertainties discussed above. The star posi- 
tions are compared graphically, one cluster at a time, with evolu- 
tion tracks for stellar masses in the range 1.7 to 9.0 M/Mq, and 
with the range of isochrones allowed by the adopted cluster or 
association age +0.1 or 0.2 dex, using models computed by the 
Padova group (see Fagotto et al. |1994| and references therein) for 
Z = 0.008. 0.02 and 0.05. (These data are available in machine- 
readable form from CDS.) This comparison allows us to test the 
hypothesis that each star is actually a cluster member, and if it 
is, to determine its mass. 

In almost all cases, the best overall agreement of individual 
star positions in the HR diagram with computed isochrones is 
found for Z - 0.02. This agreement persists even when the mag- 
netic Ap stars are near the end of their main sequence evolution. 
This Z value is somewhat higher than the currently accepted best 
estimate for the Sun, which is about Z = 0.012 (Asplund et al. 
2005). It is not clear whether our result indicates that Z w 0.02 is 
the most appropriate mean value for young stellar groups in the 
solar neighborhood, or whether this reflects a residual small sys- 
tematic error in our transformation to the theoretical HR diagram 
(for example, modest underestimation of Teff values as discussed 
above). Note that the significantly better fit of the Z = 0.02 mod- 
els to our Teff, L/Lq data compared to the fit with Z - 0.008 
models indicates that Padova models with Z = 0.012 would 
not fit our observational data as well as the Z = 0.02 models 
do. However, we cannot really test how well our data would fit 
isochrones computed with the new Asplund et al abundances, as 
these abundances have diff'erent ratios of light (CNO) elements 
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to iron peak elements than are assumed in any available grid of 
evolutionary models. 

We have tested the correctness of our computations, and as- 
sessed whether our assumed uncertainties are realistic, by look- 
ing at the results of such comparisons for nearby clusters and 
associations of well-determined membership, precisely known 
distance and moderate or small reddening. One such comparison 
is shown in Figure |2] for the Sco OB2 association, using tracks 
for Z - 0.02. This figure confirms the general correctness of the 
transformations used in going from photometric data to the the- 
oretical HR diagram, and shows that our assumed uncertainties 
give a reasonable description of the scatter about the Z = 0.02 
isochrones. 

In a few individual cases, typically in relatively distant or 
heavily reddened clusters or associations, we find significant dis- 
agreement between the deduced position of a star and the rele- 
vant isochrone which strongly suggest that the star should be 
disqualified from cluster membership on photometric grounds. 
There are eight such cases in Ori OB 1 and five in various clus- 
ters. All have been noted with "?" in the "phot" column of 
Table |2] and are omitted from Table [3] The other stars whose 
membership we have tested in this way, and whose positions are 
consistent with the isochrones to within better than 2cr are noted 
with "y" or "p" in this column. (We did not test stars whose clus- 
ter membership or Ap nature was rejected on other grounds.) 

Placing the observed stars on the appropriate range of cluster 
isochrones allows us to provide a mass estimate for each. These 
masses have been derived using the Z = 0.02 evolution tracks. 
Uncertainties for these values, also reported in Table[3] arise pri- 
marily from the size of the error oval of each stellar position 
due to uncertainties in Teff and L/Lq. A small contribution to 
the mass uncertainties is made by the uncertainty in the Z value 
appropriate for each cluster For stars whose Icr error ovals in- 
tersect the band of isochrones defined by the (inexact) age of the 
cluster, we find that the uncertainty in mass is typically about 
5 % of the actual mass value. 

Our method of mass determination may be clarified using 
Figure 12] Consider the star between the M/Mq = 6 and 7 evolu- 
tion tracks (HD 125823). It may be seen directly from the figure 
that with this set of tracks and the indicated uncertainties, the 
mass is between about 6.1 and 6.7 Mq, or M/Mq - 6.4 + 0.3. 
If the error oval of the star is a little above or below the pair of 
isochrones for the age range of the cluster, the mass uncertainty 
is increased from about 5% to about 7%. 

As noted in Sect. 3.1, there are several clusters for which 
we have determined a more precise age by forcing the cluster 
isochrone to pass through the error oval of an Ap cluster member 
near the TAMS. The method is very similar to that illustrated by 
Kochukhov & Bagnulo ( |2006| ) in their Figure 3. In these cases, 
we have included in the estimate of cluster age uncertainty a 
contribution from the uncertainty of the appropriate Z value to 
use for the isochrones. Even including this effect, in a few cases 
uncertainties in log t as small as +0.05 - 0.08 dex are obtained. 

We then use the mass and its uncertainty to derive the ex- 
pected main sequence lifetime of each star, again using the 
Z = 0.02 models. Dividing the adopted cluster age by this life- 
time, we obtain an estimate of the fraction of the main sequence 
lifetime already elapsed (the fractional age r) for each star 

Several sources of error contribute to the uncertainty in this 
quantity. Mass uncertainty of about 5 % leads to main sequence 
lifetime uncertainty of about 12% for evolution tracks for a 
given Z. The uncertainty in the actual value of Z for each clus- 
ter introduces an additional lifetime uncertainty of up to about 
7 %, leading to a total uncertainty in main sequence lifetime of 
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Fig. 3. Upper panel: Distributions of the 1 cr confidence inter- 
val of fractional age t for our sample of open cluster stars. For 
comparison, we show also the same kind of distribution for the 
sample of field stars analysed by Kochukhov & Bagnulo ( I2006I I. 
Note that Kochukhov & Bagnulo assume considerably smaller 
errors in T^ff than we do, leading to smaller uncertainties in r, as 
discussed in the text. 



about +15 %. Uncertainty in the ages of the clusters and asso- 
ciations has been estimated in most cases to be approximately 
±0.1 - 0.2 dex, or about +25 - 50%. In the clusters with age 
uncertainties in this range, cluster age uncertainty dominates the 
total error of t. In those clusters with one or more stars near 
the TAMS, for which we can obtain more precise age estimates, 
the contributions from cluster age uncertainty and from main se- 
quence lifetime uncertainty can be comparable. 

We treat the uncertainties in the stellar lifetimes and the clus- 
ter ages as independent. For stars which are young relative to 
their main sequence lifetimes (t less than about 0.3) the accu- 
racy of the fractional age is substantially better than is possible at 
present by placing field stars in the HR diagram, even for clusters 
with rather poorly known ages. For stars of fractional age below 
about 0.1, the improvement compared to the accuracy that can 
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be obtained for field stars is an order of magnitude or more. For 
T larger than about 0.5 (these are the situations in which signif- 
icant improvement in cluster age precision, by isochrone fitting 
of individual evolved stars, is possible), the uncertainty in frac- 
tional age of a cluster Ap (for given error bars on Teff and L/Lq, 
and comparable treatment of the Z uncertainty) is similar to 
that for nearby field stars, since effectively in this circumstance 
the cluster age is most precisely determined in the same way that 
the age of a field star is determined. 

The age confidence limits which characterise our sample are 
shown in Fig. [3] together with those obtained in the study based 
on field stars by Kochukhov & Bagnulo (2006). It is important to 
keep in mind when comparing these two distributions two sub- 
stantial differences between their analysis and ours. First, they 
have assumed uncertainties for T^ff that are typically about 0.012 
to 0.014 dex, compared to our uncertainty of 0.020 dex (they do 
not include possible effects due to the overall uncertainty of the 
TefF scale for Ap stars, which we consider an important probable 
source of uncertainty). Increasing their Teff uncertainty by 50% 
would increase the age uncertainty for evolved stars by a sim- 
ilar amount. Secondly, they have not included any contribution 
from uncertainty in the bulk Z values of the stars of their sample 
(and thus uncertainty in which evolution tracks to use for age de- 
termination), an effect which we have found in Paper I to make 
an important contribution to the error budget of ages for field 
stars. If Kochukhov & Bagnulo had used the same uncertainty 
estimates as we have, their 1 cr confidence limits would be very 
roughly a factor of two larger than those shown in Fig. [3] for all 
values of t. 

In spite of these important differences, it is easily seen that 
the precision of our fractional ages is considerably better than 
that of Kochukhov & Bagnulo for fractional ages below about 
0.5. For fractional ages below about 0.2, the improvement is 
about an order of magnitude, or more. Above t x 0.5, the fact 
that our fractional ages are less precise than theirs is largely due 
to the two differences in analysis noted above, since stellar main 
sequence lifetimes and cluster ages for the most evolved stars 
are determined in the same way as field star lifetimes and ages. 

The stars of the final subsample are listed in Table[3] together 
with the data derived for each star as discussed above. The first 
three columns of this table characterise the clusters whose stars 
are in our subsample, giving the cluster or association name, our 
adopted (logarithmic) cluster age log t with uncertainty, and the 
true distance modulus DM of the cluster. The remaining columns 
list the stars which we consider to be members with log Teff , 
logiL/Lo), M/Mq and uncertainty, fractional age r and uncer- 
tainty, and the RMS estimate of the longitudinal field Bims- 



5. Characteristics of our sample 

It is of interest to note some distinctive properties of the data set 
of Table [3] concerning age and mass distributions. 

5.1. Absolute and fractional age distribution 

Because open clusters and especially associations tend to be dis- 
rupted by galactic tidal acceleration on a time scale short com- 
pared to the main sequence lifetime of an A star, we have a pre- 
dominance of relatively young stars (fractional ages less than 
0.5) in our survey, although the sample does include a signif- 
icant number of more evolved main sequence stars (especially 
at higher masses where the lifetimes are shorter). In fact, about 
half the stars in the sample have fractional ages of less than 0.20. 



However, our data set is fairly uniformly distributed in (log) ab- 
solute age between log f ^ 6.5 and 8.5. The distribution of abso- 
lute ages and of fractional ages in our sample are shown in the 
first and second boxes on the left of Fig. H] respectively. 

5.2. Mass distribution 

Because the distribution of masses in a cluster or association 
follows the initial mass function (IMF), at least up to the cluster 
turn-off, rather than the field mass function, which is more heav- 
ily weighted towards low-mass stars of long main sequence life- 
times, our sample has a relatively large number of hot, high-mass 
magnetic Bp stars. The median mass in Table[3]is about 3.5 Mq, 
and the median value of Tes is about 12700 K. These are stars 
that are not particularly common among the Ap-Bp stars near 
the Sun, and may well prove to be rather interesting objects to 
study individually. The distribution in mass of stars in our sub- 
sample is shown in the third box in Fig.|4] Notice that we have 
only three stars (CPD -32 13119, CPD -20 1640, and HD 66318) 
that have M/Mq < 2.1, and of these three, only HD 66318 has 
a detected (but huge) magnetic field (cf Bagnulo, Landstreet, Lo 
Curto et al. 120031 ). Significant numbers of Ap stars with detected 
magnetic fields begin to be found in our sample only above about 

2.3 Mq. 

5.3. Magnetic field distribution 

The distribution of magnetic field strengths found for probable 
cluster Ap stars in our survey (including Ap stars in which no 
field has yet been detected) is shown in the right box in Fig. |4] 
This distribution is qualitatively similar to what is found for sam- 
ples of field Ap stars (e.g. Bohlender & Landstreet 1990), but 
scaled up to somewhat larger fields than usual. The bulk of the 
stars (about 2/3 of the present sample) have RMS fields of less 
than 1 kG, and the median field is about 450 G. This may be 
compared to the distribution of Bohlender & Landstreet; for a 
sample of 12 bright Ap field stars (also including some in which 
fields had not yet been detected) they found none with B„ns 
above 1 kG, and a median field of about 300 G. The largest dif- 
ference between this sample and that of Bohlender & Landstreet 
is the presence of a substantial tail of high-field stars; almost a 
quarter of the present sample have RMS fields of 2 kG or more. 
The differences between these two samples is probably due to 
the much larger fraction of young Ap stars in the present sample 
(see the discussion below). 



6. Discussion 

We finally consider what general conclusions may be drawn 
from our data about the evolution of magnetic fields up to and 
through the main sequence phase of stellar evolution. With a 
sample of more than 80 stars which are probable or definite clus- 
ter or association members, probable or definite Ap-Bp stars, and 
for which at least one magnetic measurement with cr of the or- 
der of 10^ G is available, we are now in a very good position 
to carry out a preliminary statistical assessment of the evolution 
of Ap-Bp star magnetic fields based on well-determined stellar 
ages. 

The fundamental results of our survey are displayed in the 
eight panels of Fig. |5] where we plot, using logarithmic scales 
on all axes, the field strength Bi^s as a function of stellar age 
(left column of the Figure), and as a function of stellar mass 
(right column). In both columns, the top panel shows all stars of 
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Fig. 4. Distributions in our sample of age, fractional age, mass, and Brm,,. In each box, the shaded histogram shows the distribution 
of those Ap stars of Table [3] for which a B„ns ^ 0.25 kG has been measured. The thicker lines show the distribution of the entire 
sample. 



Table 1. Slopes of rej 


^ression lines in Fig.|5] 




Sample 


v-axis 


X-axis 


slope 


All Masses 


log firms 


log/ 
log? 
log? 


-0.29 +0.07 
0.02 ±0.04 
-0.31 ±0.06 


M < 3.0Me 


l0g(i?2 firms) 

log/?2 

log firms 


logf 
logf 
logf 


0.07 ±0.15 
0.23 ±0.06 
-0.16 ±0.14 


3.0Me <M< A.OMo 


l0g(/;2 B,,,,) 
log firms 


logf 
logf 
logf 


-0.22 ±0.14 
0.19 ±0.06 
-0.42 ±0.14 


M > 4.OM0 


l0g(i?^ firms) 
l0g/?2 

log firms 


logf 
logf 
logf 


-0.27 ±0.12 
0.21 ±0.07 
-0.47 ±0.11 


All Ages 


l0g{R- fir„„) 

logR^ 

log firms 


log(M/Mo) 
log(M/Mo) 
log(M/Mo) 


1.44 ±0.26 
0.92 ±0.14 
0.52 ±0.29 


log t < 7.0 


logR^ 

log firms 


log(M/Mo) 
log(M/Mo) 
log(M/Mo) 


0.84 ±0.57 
1.12 ±0.21 
-0.29 ±0.61 


7.0<logf<8.0 


log(R^ B,^,) 
logR^ 

log firms 


log(M/Mo) 
log(M/Mo) 
log(M/Mo) 


1.23 ±0.52 
1.32 ±0.14 
-0.09 ±0.51 


log t > 8.0 


l0g(R'- firms) 

logR^ 

log firms 


log(M/Mo) 
log(M/Mo) 
log(M/Mo) 


0.70 ±0.61 
1.83 ±0.38 
-1.13 ±0.57 



the sample (except Ori C), while below we see the same plot 
for three separate mass ranges (left column) or age ranges (right 
column), with parameter ranges shown in each panel. 

The use of logarithmic scales would introduce a significant 
bias into fits to these data, in that an RMS field which happens to 
have a value of only a few G is substantially lower than the point 
for a star with an RMS field of, say 150 G, although with the un- 
certainties of our data, there is no statistical difference between 
the two values. In order to avoid having small RMS field values 



(essentially at present, non-detections) have an undue influence 
on fits, we have artificially offset all firms values by H-200 G. 

6.1. Does magnetic fields strength change with with time? 

The fundamental effect we are looking for with this study is 
to discover whether typical RMS fields of magnetic Ap stars 
change with stellar age, and if so, how they change. The top left 
panel of Fig. |5] appears to provide a first tentative answer to this 
question. The variation of log Bmis as a function of log t has been 
fit with a linear function. This fit to the (modified) data, shown as 
a line through the data points, was found to have a slope which 
is non-zero (negative) at about the 5cr level (see Table [T]i. It ap- 
pears that the typical field strength decreases markedly with stel- 
lar age, from a value of more than lOOOG for stars having ages 
of about 3 10'' yr to a value of the order of 200 G at 6 10** yr. 

However, from Table [3] it is clear that most of the stars with 
very young ages (say, less than 3 10^ yr) are relatively mas- 
sive (masses around 4 or 5 Mq), while the stars with ages of 
more than 10^ have smaller masses (below about 4 Mq). Another 
symptom of this systematic change in mass with age is seen 
in the variation of radii (in fact, we have plotted the value of 
(R/Rq)^) of the stars: the radius of a single star increases with 
age, but the radii of the sample show no significant variation 
with age, because the increase in radius of individual stars is es- 
sentially compensated by the fact that older stars have smaller 
masses and radii. 

Thus we have considered three smaller sub-samples having 
limited mass ranges: M < 3 Mq (shown in the second panel from 
the top on the left), 3 Mq < M <4 Mq (third panel from the top), 
and M > 4Mq (left bottom panel). Note that the subsample of 
stars defined by M < 3 Mq is in fact mainly populated with stars 
between 2.3 and 3 Mq, and the subsample defined by M > 4, Mq 
in fact contains mainly stars between 4 and 5 Mq. The selected 
mass ranges are narrow enough that the inhomogeneity of mass 
is not a major factor 

It is seen that for the two samples with larger masses (3 - 
4 Mq and > 4 Mq) the Z?r,ns decreases with time with a slope that 
is non-zero at the 3 cr level. For the 3-4 Mq sample, more than 
half of the stars younger than 410^ yr have RMS fields greater 
than 1 kG, while all but one of the stars in this range older than 
this age have fields below 1 kG. Similarly, in the most massive 



J. D. Landstreet, S. Bagnulo, V. Andretta, et al.: The evolution of magnetic fields in Ap stars 



11 



10^ 
104 

103 
102 

101 
104 

S 103 



102 

S 101 

6 104 



103 
p 102 

^ 101 

104 
103 
102 
101 




All Masses 







c 









o 


o 




^ g 


o 





M < 3.0 M 




= 



108 




I M > 4.0 M 



All Ages • T 



10= 




Age (years 



4 5 6 7 

Mass (M ) 



Fig. 5. Plots of Brms (stars), (R/Rq)^ (open circles), and B„ns(^/^Q)^ (filled circles) are shown as functions of stellar age (left 
column) and mass (right column) for the subsample of Table [3] In each column the top panel shows all stars in the sample (except 
for 6^ Ori C). Lower panels show various cuts through the subsample as indicated on each panel. As discussed in the text, 200 G 
has been added to each value of B„ns, in order to avoid having correlations made artificially steeper by very small values of firms- 



sample, a majority of the stars with ages below 3 10^ yr have 
fields above 1 kG while none of the older stars (with a range of 
ages up to more than 10** yr) have fields above 1 kG. 

In contrast, the lowest mass sample (mainly populated with 
stars between 2.2 and 3 Mq) does not show a significant decrease 
of RMS field with age, and in fact there are about as many fields 
above 1 kG among the older (age > 10^ yr) stars as among the 
younger stars of this mass range. Note that this group has no rep- 
resentatives among the youngest clusters; we have not found any 



magnetic Ap stars with M < 3Mq in stellar groups of ages less 
than 10^ yr, although the Sco OB2 association is close enough 
for such stars to be readily identified. 

One possible simple scenario for the observed decline of 
firms with age among stars of M > 3Mq would be that such stars 
(approximately) conserve magnetic flux as they evolve through 
the main sequence phase. During the evolution from ZAMS to 
TAMS, the radius of a star increases by a factor of roughly 3. 
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We have tested our samples to see if flux Ob ~ Bi-^^iR/Rof- 
is statistically conserved as we look at stars of various ages. In 
the various panels of Fig. |5] we plot both (R/Rq)^ and B^^sR^- 
In the left column of panels, the individual mass bands all show 
increasing values of (R/Rq)^ with increasing age at the 3cr level, 
although the rise is not as great as expected, probably due to the 
increasing contribution of lower mass stars in each sample as we 
go to greater ages. In the two panels at the lower left, where we 
find strongly significant decline in Bi^s, the value of flux Brms^^ 
appears to decline somewhat, but the decline is only present at 
about the 2cr level. We may regard this as consistent either with 
some real decline in magnetic flux, or with conservation of mag- 
netic flux. 

To summarise the results of this analysis, we find that mag- 
netic Ap stars with masses above 3 Mq have fields which decline 
substantially over an age of the order of 2 - 3 10^ yr For these 
stars, the magnetic flux may decline with time, or remain roughly 
constant. The RMS fields of less massive stars change little over 
time-scales approaching lO'' yr, and there is no strong evidence 
that the magnetic fluxes change significantly over this time scale. 

The evolution with time of fields in main sequence Ap stars 
has also been discussed by Kochukhov & Bagnulo (2006 ). They 
find that the value of B„ns declines with increasing age for all 
their subsamples of Ap stars except for those with M/Mq < 2, 
in contrast to our result that B^^^ decreases with age except for 
stars of M/Mq < 3. Thus the two studies agree that there is a sig- 
nificant decrease of Bi^s with time among Ap stars at least above 
3Mq, although they disagree on precisely where this behaviour 
changes over to more nearly constant fields. 

However, Kochukhov & Bagnulo (2006) find a significant 
increase in magnetic flux with age for stars of M/Mq < 3, while 
we find, if anything, a modest decrease of flux with time. At this 
point, it is probably fair to say that approximate magnetic flux 
conservation during the main sequence phase is not ruled out. It 
will be very useful to obtain more precise values of Bims for stars 
in the cluster sample (by getting multiple observations of each 
star, and expanding the sample of stars measured) to resolve the 
differences between these two studies. 

6.1 .1 . Do magnetic fields appear late in main sequence 
evolution? 

So far we have discussed a possible decline of the magnetic field, 
but we have not commented on whether the magnetic field is 
present at the very moment when a star reaches the ZAMS, or 
whether the field appears at the stellar surface at some later stage 
in its main sequence evolution. This problem has been previ- 
ously discussed, among others, by Hubrig et al. (2000), Bagnulo 
et al (I2003l l. Pohnl et al. (2005), and Kochukhov & Bagnulo 
(I2OO6I 1. 

Based on the analysis of field stars, Hubrig et al. ( I2000I I 
have proposed that, in stars with M ^ 3 Mq, a magnetic field 
appears at the surface of an Ap stars only after about 30 % of 
its main sequence liftime has elapsed. This conclusion was later 
reaffirmed by Hubrig et al. (|2005al|2005bl l. However, this result 
was contradicted by the discovery of a field in the young 2.1 Mq 
star HD 66318 (Bagnulo et al. |2TO3] Furthermore, both Pohnl et 
al. ([2005 ) and Kochukhov & Bagnulo ( |2006l l were unable to re- 
produce the results of Hubrig et al. of late emergence of surface 
fields, but both studies found a remarkable shortage of young 
magnetic stars with M ^ 2 Mq. 

In the present sample we have a total of six stars with masses 
definitely below 3Mq, fractional ages below 0.30, and unam- 
biguously detected magnetic fields. These stars are HD 65712 



M < 3M,. 
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Fig. 6. The position of the stars with M < 3 M© in the diagram 
of firms vs. fractional age. 



(Paper I), HD 66318 (Bagnulo et al l2003l l. HD 92385 (Paper 
I), HD 1 12381 (Bohlender et alT99j|, HD 1 19419 (Thompson, 
Brown & Landstreet 1987 ), and HD 318100 (Paper I). One addi- 
tional star, HD 74169 (Paper I), satisfies all criteria but only has 
one measurement with field detection at the few sigma level, and 
still ought to be confirmed. Figure |6] shows the position of all of 
the stars of our sample with M < 3 Mq in the diagram of Z?i,ns 
vs. fractional age. Our study thus definitely contradicts that of 
Hubrig et al. concerning late field emergence. 

We have investigated the reasons for the discrepancies be- 
tween our conclusions and those of Hubrig et al. (I2000I I and we 
have identified several points that may have led them to obtain 
incorrect results. 

i) Hubrig et al. (120001 ) have used the bolometric corrections 
of Schmidt-Kaler, which are too large even compared to the most 
accurate values for normal stars, and lead to an overestimate 
of bolometric magnitude relative to the BC's adopted here (see 
Fig. [TJ by about 0.25 mag (0.1 dex in log(L/L0)). Using more 
appropriate BC's would move their observed stars downwards, 
towards the ZAMS, and younger ages. 

ii) As discussed above, the whole effective temperature scale for 
magnetic Ap stars may be systematically a few percent low. This 
possibility has little effect on ages derived for cluster stars, but, 
if correct, would move the field stars towards the ZAMS as well. 
Hi) Hubrig et al. (2000) applied the Lutz-Kelker correction to 
their data, moving the stars upwards in the HR diagram, away 
from the ZAMS. As discussed by Stgpien ( 120041 ) and Kochukhov 
& Bagnulo (2006), this correction should not have been used. 
Removing it also moves the Hubrig et al. sample closer to the 
ZAMS. 

iv) Finally, in Hubrig et al. (120001 ) there are a two significant er- 
rors: the Teif value of HD 1 19419 is about 15 % too low, and that 
of HD 147010 is more than 50 % too small. Correcting these er- 
rors would move HD 147010 to a mass slightly above 3Mq but 
to a position close to the ZAMS. HD 119419 would have re- 
mained well below 3 Mq and moved close to the ZAMS, provid- 
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ing Hubrig et al. with a counter-example to their hypothesis in 
their own data. 

The changes to the data of Hubrig et al. (I2000I I resulting from 
the points discussed above appear to be enough to make the ef- 
fect they report disappear from their data. This result, together 
with the results of Bag nulo et al. (I2003I I. Pohnl et al. (2005), 
Kochukhov & Bagnulo ('2006), and our own survey, lead us to 
reject the hypothesis that stars with M < 3 Mq becomes mag- 
netic only after they have spent a significant fraction of their life 
in the main sequence. 

In our survey we note a relatively small, but non-zero, num- 
ber of stars of mass below about 2Mq (T^ff < 8500 K). Only 
three of the 81 stars in Table [3] have masses in the range 1.7 < 
M/Mq < 2.1. To some extent this may reflect the absolute rar- 
ity of such stars; in the field, the Ap stars near FO make up only 
about 1 % of all stars, compared to about 10 % near AO (Wolff 
I19681 I. Still, a fairly large number of such stars are now known 
in the field, both as a result of searches for roAp stars (Kurtz & 
Martinez 119931 ) and from searches for magnetic stars showing 
visible Zeeman splitting, which seems to be commoner among 
low-mass Ap stars (Hubrig et al. 2000). A number of such low- 
mass Ap stars are also found in the survey of Kochukhov & 
Bagnulo (2006 ); the rarity of low-mass Ap stars in our sample 
seems to be anomalous. It may result primarily from selection 
effects. Such stars will be among the faintest Ap stars in most 
of the clusters we have studied, and may simply not have been 
selected as Ap stars in identification studies, either because they 
were near the threshold of the study, or because the method used 
(such as the Aa method) is not very sensitive to the coolest Ap 
stars. 

Note that, in contrast to the results of Pohnl et al. ( |2005| l and 
Kochukhov & Bagnulo (120061) . both of whom found a deficiency 
of young Ap stars with masses below 2Mq, all three of the low- 
mass stars in our survey are fairly close to the ZAMS; the largest 
fractional age among these stars is 0.16. This effect probably 
results from the fact that the main sequence lifetime of such stars 
(10^ yr or more) is considerably longer than the ages of most 
clusters in our sample. Another distinctive feature of our sample 
of particularly low-mass Ap stars is that no field has so far been 
definitely detected in either of the two stars with masses below 
2Mq. The lowest mass cluster Ap star for which a field has been 
definitely detected is HD 66318 (Bagnulo et al. |2003] ). This is 
quite different from the situation in the field, where a number of 
Ap stars with detected fields have masses that extend down to 
about I.5M0 (cf Kochukhov & Bagnulo). 

From the present study (in particular from the fractional ages 
reported in Table [3j it is clear that magnetic fields are present 
at essentially the zero age main sequence phase for all stellar 
masses between about 2 and 5Mq; stars with fractional ages be- 
low about 0.05 occur throughout this range. It seems very rea- 
sonable to suppose that these fields are already present late in the 
pre-main sequence phase, and this has recently been confirmed 
by the discovery of ordered global fields in several pre-main se- 
quence Her big Ae Be stars (Donati et al. |1997t Wade et al. i2005i 
Catala et al. |IU06l ). 

6.2. Do magnetic fields depend on mass? 

A further issue that we may address with our data is the question 
of whether there is a variation with mass of typical RMS field 
strength, at a given stellar age. Thompson, Brown & Landstreet 
(|1987 ) have argued that they detected such an effect, with stars 
more massive than roughly 6 or 7 Mq having fields about a fac- 
tor of two larger than lower mass stars. Kochukhov & Bagnulo 



( 120061 ) find that the field strength Brma shows no clear variation 
with mass, but argue that the magnetic flux definitely increases 
with mass. 

A plot of the RMS fields of all the stars in our sample as 
a function of stellar mass (top right panel in Fig. |5]l appears to 
show a trend of increasing field strength with mass. However, 
again we must recall that this sample has considerable change in 
typical mass as a function of age. To reduce the influence of this 
effect, the whole sample is subdivided into three samples of more 
Hmited age range: log t < 7.0, 7.0 < log f < 8.0, and log t > 8.0. 
As noted above, the youngest sample contains only stars more 
massive than about 3 Mq, and the oldest sample is restricted to 
stars less massive than 4 Mq. None of these smaller samples of 
restricted age shows a clearly significant trend of RMS field with 
mass. In particular, we have six stars in the youngest sample with 
masses above the value at which Thompson et al. ( 119871 ) thought 
they had found a jump in field strength, but we see no significant 
change at this mass. 

From our data, it is not clear if magnetic flux is typically a 
function of stellar mass. The age range between log f = 7 and 
8 appears to show a significant increase of magnetic flux with 
mass, but this does not seem to be present (over more restricted 
mass ranges) for younger or older samples. 

7. Conclusions 

We have assembled a sample of nearly 160 stars for which mag- 
netic field measurements are available (primarily from Paper I) 
and which are possible members of clusters or associations, and 
possible magnetic Ap or Bp stars. After detailed examination 
of astrometric, photometric, and spectroscopic data available for 
these stars, 81 remain as probable cluster members and probable 
magnetic Ap stars. These stars are listed in Table [3] 

For these stars, we have obtained cluster/association dis- 
tances and ages (typically accurate to +0.05 - 0.2 dex) from 
the literature and from our own isochrone fits to stars near the 
TAMS. We have determined effective temperatures (typically 
accurate to +0.02 dex) for almost all of these stars from both 
Stromgren and Geneva photometry. We have derived new bolo- 
metric corrections for Ap stars, and obtained luminosities of in- 
dividual stars with typical uncertainties of ±0. 1 dex. 

By placing the stars of this subsample in the HR diagrams 
of their clusters, we have determined their masses with an ac- 
curacy of about ±5 %, and improved the precision of the ages 
of several clusters. The mass values allow us to obtain main se- 
quence lifetimes from evolutionary calculations, with an overall 
uncertainty (including the uncertainty in bulk chemical composi- 
tion) of about ±15 %. From main sequence lifetimes and cluster 
ages, we deduce fractional ages t (the fraction of the main se- 
quence lifetime of each star already elapsed). The uncertainty 
in the fractional age is almost always less than half of the actual 
fractional age, and thus for stars near the beginning of their main 
sequence lives, the fractional ages are much more accurate than 
ages than can be derived for field stars by placing them in the 
HR diagram (see the penultimate columns of Table[3ll. 

From the stellar characteristics derived here we have as- 
sembled for the first time a substantial sample of magnetic Ap 
stars for which the absolute and fractional ages are known with 
really useful precision. These accurate masses and fractional 
ages allow us to definitively test the hypothesis of Hubrig et 
al. ( 12000 1) that in Ap stars of masses M < 3Mq, the magnetic 
fields first become visible at the surface after about 30 % of the 
main sequence lifetime has elapsed. We have identified six clear 
counter-examples to this hypothesis in our sample. As a result. 
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we are able to definitively rule out the hypothesis of late field 
emergence in low-mass Ap stars. 

Although the field measurements of the stars of our sample 
are still rather incomplete, we can use the available data to con- 
strain the typical evolution of magnetic field strength with stellar 
age. We find that fields are definitely present essentially at the 
ZAMS (fractional ages of 0.05 or less) for Ap stars ranging be- 
tween 2 and 5 Mq, and for fractional ages of less than about 0. 10 
for masses up to 9 Mq. 

The evolution of magnetic fields with time suggested by our 
data reveals a quite unexpected pattern. For stars of M > 3Mq we 
find strong evidence that the field strength declines by a factor 
of a few on a timescale of about 2 - 3 10^ yr, in agreement with 
the conclusion of Kochukhov & Bagnulo ( |2006l l. It is not clear 
for these stars whether the total magnetic flux remains approx- 
imately constant, causing B^ms to decrease as the stellar radius 
expands, or whether this flux actually declines slowly. 

In contrast, for stars of M < 3Mq, there is no convincing 
evidence in our data of field strength decrease, nor of mag- 
netic flux decrease, even on a time scale of several times 10** 
yr. This conclusion, very different from the result for more mas- 
sive Ap stars, represents a significant puzzle; perhaps this result 
reveals some important structural difference between the lower 
and higher mass Ap stars. 

Finally, we are able to examine the question of whether there 
is a significant variation in typical magnetic field strength as a 
function of mass within a given age band. We do not find any 
strong evidence for such a variation, in contrast to Thompson et 
al. (fl987), who argued that typical field strength increases rather 
abruptly at about 6-7 Mq. 

These results indicate clearly that studies of magnetic stars 
in clusters can yield much useful information about the nature 
and evolution of their fields. Further observations, both to detect 
fields in more cluster Ap stars (thereby increasing the size and 
diversity of the sample), and to improve the accuracy of the RMS 
fields in stars already detected, are already in progress. 
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Table 2. Tests of cluster membership and chemical peculiarity. The first two columns give the cluster or association to which the star may belong, 
and a common designation of the star. Columns 3, 4 and 5 indicate whether the measured parallax, proper motions, and position in the HR 
diagram are consistent with membership, and column 6 gives our conclusion about cluster membership (y = yes, p = probable, ? = questionable, 
n = no, blank = no data; letters in parentheses refer to references at end of table). Columns 7 through 10 summarize evidence that each star is an 
Ap star according to its spectral classification, to the value of Aa and/or Z, to the presence of periodic variability, and to whether a magnetic field 
is actually detected (letters in parentheses in the "mag fid" column refer to references at the end of the table). Column 1 1 reports our conclusion as 
to whether the star is indeed an Ap star. 
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Notes on individual stars (* in "member?" or "Ap?" column): 

HD 19805 Composition normal, no field, not Ap star (private communication, J. Silvester) 

HD 23387A Composition normal, no field, not Ap star (private communication, J. Silvester) 

HD 37022 Tycho-2 jis inconsistent with Hipparcos value and with membership of this Trapezium star in Ori OB Id. 

Accept Hipparcos value. 

CpD -20 1640 This star is Cox 40. Simbad identification and coordinates of NGC 2287 40 are incorrect. 

HD 56343 Parallax and proper motions agree with cluster means, but Baumgardt et al 00 regard this star as 

non-member We accept it as member 
HD 62376 Canier et al. <1999t and Baumgardt et al. 420001 accept star as probable inember Because of 5cr discrepancy in 

fig we consider it a non-member 
HD 65712 Star seems to be outside field considered by Dias et al 2001. 

CpD -60 944B Dias et al 2001 consider this star non-member, but both proper motions differ from cluster mean by less than 

2<T. Accept as probable member 
HD 65949 Both Hipparcos and Tycho-2 proper motions have large uncertainties, but are consistent with probable 

membership. Disregard small inembership probability assigned by Dias et al 2001. 
HD 74195 Hipparcos and Tycho-2 jU^'s inconsistent. Accept Tycho-2 value which is consistent with membership. 

HD 136347 Tycho-2 /j's both completely different from Hipparcos values. Accept Hipparcos values and membership. 
HD 145501 Membership assessment for v Sco CD based on tt and /j's for v Sco AB. 

HD 190290 Cluster mean motions unknown for Melotte 227. 



References for "member?" column: 


(B) 


Baumgardt et al. ('2000) 


(Br) 


Brown et al. ( 1999i 


(C) 


Carrier et al. 1 1999j 


(D) 


Dias etal. (2001 1 


(D2) 


Dias et al. 42002) 


(K) 


Kharchenko et al. (2005t 


(M) 


Marschall etal. 41982) 


(R) 


Robichon et al. 4 19991 


(S) 


Sabogal-Martinez et al. 4200U 


(T) 


Tian et al. 41996) 


(Z) 


deZeeuwetal. 41999) 



References for "mag fld" column: 



(F) 


Bagnulo et al. (2006 1, FORSl data froin Paper 1 


(a) 


Bychkov et al. t2003( 


(b) 


Brown etal. 419851 


(c) 


Borra 41981) 


(d) 


Bohlender et al. 41987) 


(e) 


Wade et al 42006) 


(g) 


Thompson & Landstreet 41985) 


(h) 


Bohlender et al. 41993) 


(i) 


Preston etal. 41969) 


G) 


Shorlin et al. 42002) 


(k) 


Bon-aetal. 41983) 


(1) 


Landstreet 41990) 


(m) 


Landstreet et al. ( 19^ 


(n) 


Thompson et al. (1987) 


(o) 


Mathvsetal. U997) 


(P) 


Hubrig et al. )2004| 


(q) 


Babcock ( 1958 1 


(r) 


unpublished Musicos data, courtesy J. Silvester 


(s) 


Kochukhov & Baenulo J2006l 



J. D. Landstreet, S. Bagnulo, V. Andretta, et al.: The evolution of magnetic fields in Ap stars, Online Material p 5 
Table 3. Physical properties of stars that are probable open cluster members, probable Ap stars, and have magnetic field measurements. 
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0.30 ± 0.09 


610 








HD 146555 


4.107 


2.19 


3.45 ± 0.20 


0.33 ±0.10 


270 


NGC 6178 


7.15 ±0.22 


10.03 


HD 149257 


4.406 


3.82 


10.30 ± 0.60 


0.61 ±0.33 


160 








HD 149277 


4.348 


3.71 


8.75 ± 0.40 


0.45 ± 0.24 


2435 


NGC 6193 


6.69 ± 0.20 


10.31 


CD-48 11051 


4.301 


3.32 


7.00 ±0.30 


0.10 ±0.05 


2400 


NGC 6281 


8.45 ±0.10 


8.86 


HD 153948 


4.025 


1.86 


2.70 ±0.15 


0.62 ±0.18 


195 


NGC 6405 


7.80 ±0.15 


8.44 


HD 318107 


4.072 


1.92 


2.95 ±0.15 


0.17 ±0.07 


4820 
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Cluster 


\ogl 


true DM 


Star 


logr. 


log L/Lq 


M/Mo 


fractional age 










HD 318100 


4.025 


1.64 


2.50 ±0.10 


0.11 ±0.04 


345 








CoD -32 13119 


3.892 


1.06 


1.75 ±0.10 


0.04 ± 0.02 


35 


NGC 6475 


8.47 ±0.13 


7.24 


HD 162305 


4.004 


1.82 


2.65 ±0.15 


0.61 ± 0.22 


85 








HD 162725 


3.982 


2.36 


3.30 ± 0.20 


1.08 ±0.38 


60 


NGC 6633 


8.75 ± 0.05 


7.93 


HD 169842 


3.924 


1.62 


2.35 ±0.10 


0.85 ±0.18 


370 


IC 4725 


8.02 ± 0.08 


10.44 


BD-19 5044L 


4.107 


2.25 


3.55 ±0.15 


0.46 ±0.12 


235 








HD 170836 


4.133 


2.80 


4.55 ± 0.20 


0.84 ± 0.22 


505 








HD 170860 


4.137 


2.63 


4.25 ± 0.20 


0.71 ±0.19 


40 
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